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ABSTRACT 


Experimental  work  was  carried  out  on  the  steady  state  heat  transfer  behavior  from 
a  tube  bank  in  a  zero  mean  oscillatory  flow.  The  oscillatory  flow  across  the  tube  bank 
was  created  by  an  acoustic  field  inside  an  isolated  resonant  chamber.  The  transverse  tube 
bank  arrangement  was  represented  by  smooth  walled  cylinders  placed  parallel  to  each 
other,  with  the  plane  of  the  cylinders  normal  to  the  direction  of  fluid  oscillation,  similar 
to  the  arrangement  found  in  many  heat  exchangers.  The  spacing  between  the  cylinders 
was  varied  to  examine  the  effects  of  boundary  layer  interference  on  the  heat  transfer 
behavior  .  Heat  transfer  correlations  were  developed  in  the  form  of  Nusselt  number  as  a 
function  of  the  streaming  Reynolds  number  for  each  tube  spacing.  This  experimental 
study  is  relevant  to  the  design  of  heat  exchangers  for  thermoacoustic  engines. 
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I.  INTRODUCTION 


The  influence  of  sound  waves  in  the  transfer  of  heat  is  a  relatively  new  science.  A 
recent  application  of  acoustic  heat  transfer  is  relevant  to  the  design  of  heat  exchangers  for 
thermoacoustic  engines/refrigerators.  Unlike  the  vapor  cycle  refrigerators  in  widespread 
use  today,  thermoacoustic  refrigeration  does  not  involve  the  use  of  chloroflourocarbons 
(CFCs)  that  are  potentially  harmful  to  the  ozone  layer  in  the  atmosphere.  Additionally, 
thermoacoustic  refrigerators  do  not  contain  many  moving  parts,  which  can  substantially 
reduce  maintenance  costs.  Thermoacoustic  refrigeration  uses  sound  waves  to  remove  heat 
from  the  refrigerated  space  to  produce  the  cooling  effect.  Unlike  uniform  flow,  however, 
there  is  little  data  available  to  quantify  the  heat  transfer  exchange  with  bodies  in  an 
oscillating  flow. 

In  order  to  make  the  commercial  use  of  the  thermoacoustic  refrigerator  an 
eventual  reality,  it  must  be  economically  competitive  with  conventional  vapor  cycles.  At 
present,  thermoacoustic  refrigeration  is  less  efficient  than  conventional  cooling. 

Designing  effective  heat  exchangers  for  use  in  an  acoustic  field  is  a  crucial  step  in 
making  these  refrigerators  more  efficient.  To  accomplish  this,  experimental  data  must  be 
collected  to  correlate  the  heat  transfer  rate  from  a  heat  exchanger  to  the  acoustic  field 
characteristics  so  as  to  allow  a  more  knowledgeable  design. 

Wheately  et  al.  (1983)  and  Atchley  et  al.  (1990)  have  performed  detailed 
experiments  to  examine  the  heating/cooling  thermoacoustic  effect  produced  by  a  standing 
acoustic  wave.  Raney  et  al.  (1954),  Westervelt  (1960),  and  Nyborg  (1965)  have 
developed  theory  for  the  behavior  of  acoustic  streaming  flow  around  isolated  bodies. 
Mozurkewich  (1995),  and  Gopinath  and  Harder  (2000),  have  performed  heat  transfer 
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experiments  with  an  isolated  cylinder  in  an  oscillating  flow.  Correlations  to  relate  the 
dimensionless  parameters  important  to  heat  transfer,  namely  the  Nusselt  Number  (Nua)  as 
a  function  of  the  Streaming  Reynolds  Number  (Rs)  were  found  by  Gopinath  and  Harder 
(2000).  Common  heat  exchanger  designs,  however,  usually  consist  of  not  a  single 
isolated  cylinder,  but  several  tubes  arranged  in  various  geometries.  One  common 
geometry  is  a  bank  of  tubes  in  a  transverse  arrangement,  i.e.  placed  with  its  plane  normal 
to  the  direction  of  flow  -  it  is  this  particular  geometry  which  forms  the  subject  of  study  in 
this  thesis.  Depending  on  the  distance  between  the  tubes,  the  boundary  layers  formed  on 
the  surface  of  the  tubes  by  the  oscillating  flow  may  interfere  with  each  other  and  have  a 
deleterious  affect  on  the  heat  transfer  rate.  It  would  be  beneficial  then,  to  determine  the 
extent  of  this  interference  so  that  heat  exchanger  design  could  incorporate  this  influence. 

This  experimental  work  involves  the  correlation  of  the  Nusselt  number  to  the 
streaming  Reynolds  number  for  a  transverse  bank  of  cylindrical  tubes  subjected  to  an 
acoustic  field.  A  high  power,  standing,  resonant  wave  is  generated  in  an  isolated  chamber 
by  an  acoustic  compression  driver.  The  resulting  oscillating  cross  flow  is  used  to  remove 
heat  from  a  heated  cylinder  with  unheated  or  “dummy”  cylinders  placed  alongside  it  to 
simulate  the  various  neighboring  tubes  as  in  a  heat  exchanger  bank.  The  data  required  to 
calculate  the  Nusselt  number  and  streaming  Reynolds  number  are  gathered  using  various 
electronic  instruments.  The  uniform  spacing  between  the  cylinders  is  parametrically 
varied  so  that  the  resultant  effect  of  the  boundary  layer  interference  on  the  heat  transfer 
can  also  be  determined. 
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II.  BACKGROUND 


A.  ACOUSTIC  FIELDS 

Classical  fluid  dynamics  and  convective  heat  transfer  studies  have  for  obvious 
reasons  focused  most  attention  on  uniform  fluid  flows,  i.e.  flows  which  do  not  reverse 
direction,  let  alone  do  so  in  a  rapid  and  repetitive  manner  as  in  an  acoustic  field.  The  fluid 
is  put  into  motion  by  some  force,  such  as  a  fan,  pump,  or  atmospheric  wind.  The  moving 
fluid  then  transfers  heat  into,  or  away  from,  a  surface,  and  is  then  carried  away.  A 
common  example  is  the  air  flowing  through  the  finned  passages  of  the  radiator  of  an 
automobile.  This  situation  could  be  modeled  as  a  uniform  flow  of  cool  air  at  a  velocity 
approximately  equal  to  the  speed  of  the  vehicle,  flowing  across  finned  tubes  carrying  hot 
engine  coolant. 

Oscillating  flow,  on  the  other  hand,  is  a  situation  in  which  the  fluid  flow  direction 
changes  continuously  back  and  forth.  The  fluid  velocity  has  a  zero  mean  since  the 
amplitudes  of  oscillation  in  both  directions  are  equal,  but  their  directions  are  opposite. 
There  is  no  net  movement  of  the  fluid  over  time.  However,  there  can  still  be  heat  transfer 
in  an  oscillating  flow  because  convection  occurs  due  to  the  basic  motion  of  the  fluid, 
regardless  of  the  fluid  direction.  Heat  transfer  from  an  oscillating  flow  has  been  studied 
more  and  more  frequently  in  the  latter  half  of  the  twentieth  century,  as  reviewed  for 
instance  by  Richardson  (1967). 

In  an  acoustic  sound  field,  an  acoustic  driver  vibrates  and  causes  the  fluid 
oscillation.  If  a  resonant  field  is  isolated  inside  a  chamber,  such  as  a  plane-ended 
cylinder,  the  field  will  have  a  standing  wave  distribution,  i.e.  sinusoidally  varying 


pressure  and  velocity  distributions  that  are  90  degrees  out  of  phase  with  each  other.  A 
pressure  antinode,  or  maximum,  is  created  at  the  planar  rigid  end  termination  of  the 
chamber  where  the  velocity  is  zero  from  the  no-slip  boundary  condition.  This  situation  is 
illustrated  in  Figure  1.  Since  the  pressure  and  velocity  fields  are  90  degrees  out  of  phase 
(like  the  sine  and  cosine  functions),  the  pressure  will  ideally  be  at  an  antinode  where  the 
velocity  is  at  a  node,  and  vice  versa.  This  characteristic  has  major  implications  as  far  as 
heat  transfer  is  concerned,  as  will  be  discussed  later. 


Rigid  Plane 
Terniinrtion 


Acoustic 

Dxivtx 


Figure  1:  Resonate  Standing  Acoustic  Wave  in  an  Isolated  Chamber 


B.  THERMO  ACOUSTIC  REFRIGERATION 

A  temperature  gradient  along  a  tube  can  cause  pressure  oscillations  in  the  fluid 
along  the  tube  surface,  which  can  create  audible  sound  under  the  right  circumstances.  The 
basic  principle  behind  thermoacoustic  refrigeration  is  the  fact  that  this  process  can  also  be 
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reversed,  i.e  an  acoustic  field  under  the  right  circumstances  can  be  used  to  create  a 
thermal  effect.  Garrett  and  Hofler  (1992),  and  Swift  (1995),  describe  thermoacoustic 
refrigeration  in  detail.  The  simplified  model  in  Figure  2  depicts  a  basic  thermoacoustic 
refrigerator.  A  stack  of  thin  plates  with  low  thermal  conductivity  is  placed  at  a  specific 
location  within  an  isolated  sound  chamber.  An  acoustic  driver  at  one  end  oscillates  the 
gas  (typically  air  or  a  low  Prandtl  number  gas  mixture)  inside  the  chamber.  The  fluid 
oscillates  between  high  and  low  pressure  locations  inside  the  tube.  According  to  classical 
thermodynamics,  as  a  gas  particle  is  compressed  in  a  high-pressure  region,  its 
temperature  increases.  It  will  therefore  give  heat  off  to  the  stack  at  that  location.  As  it 
oscillates  back  to  a  low-pressure  region  and  expands,  its  temperature  decreases,  and  it 
removes  heat  from  the  stack  at  that  location.  A  heat  exchanger  at  each  end  of  the  stack  is 
where  the  stack’s  temperature  gradient  is  utilized  for  refrigeration.  Essentially,  sound  is 
the  mechanical  energy  that  serves  to  pump  heat  from  the  cold  exchanger  to  the  hot 
exchanger.  The  magnitude  of  the  temperature  gradient  along  the  stack  increases  as  the 
pressure  ratio  of  the  fluid  oscillation  increases. 


5 


Figure  2:  Thermoacoustic  Refrigerator 

C.  BOUNDARY  LAYERS 

To  utilize  the  temperature  gradient  along  the  stack,  heat  exchangers  must  be  able 
to  transfer  heat  away  from  the  low  temperature  end  of  the  stack  to  the  high  temperature 
end.  These  heat  exchangers  must  be  able  do  this  as  efficiently  as  possible  if  a  large 
temperature  gradient  is  to  be  obtained. 

In  convective  heat  transfer  with  a  uniform  flow,  a  higher  fluid  velocity  will 
generally  mean  a  higher  convective  heat  transfer  coefficient  (h).  The  same  is  true  in  an 
oscillating  fluid  flow.  It  would  seem  logical,  therefore  to  place  the  heat  exchangers  at 
pressure  nodes,  where  the  velocity  is  at  a  maximum.  This  condition  may  require  that  a 
gap  in  physical  contact  exist  between  the  stack  and  the  heat  exchanger.  It  is  therefore  all 
the  more  important  for  the  heat  exchanger  to  be  efficient. 

The  most  common  heat  exchanger  designs  utilize  several  cylinders  placed  parallel 
and  in  line,  with  their  plane  perpendicular  to  the  cross  flow  of  the  fluid.  Garrett  et  al. 
(1994)  have  provided  a  preliminary  heat  exchanger  design  based  however  on  the 
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principles  of  conventional  mean  flow  analysis.  Figure  3  shows  a  common  heat  exchanger 
tube  bank  type  arrangement. 


Figure  3:  Heat  Exchanger 

If  the  parallel  cylinders  in  a  flow  are  close  enough  to  each  other,  the  boundary 
layer  flow  formed  around  each  cylinder  could  interfere  with  the  boundary  layer  on  the 
adjacent  cylinders.  Figure  4  shows  this  situation.  For  this  boundary  layer  interference  to 
occur,  the  boundary  layer  thickness,  (6),  must  be  at  least  half  of  the  distance  between  the 
cylinders,  (Q.  Therefore, 

2  S>C 

is  the  requirement  for  boundary  layer  interference  to  occur.  In  oscillating  flows,  8 
depends  on  several  parameters,  including  the  frequency  of  fluid  oscillation,  as  will  be 
discussed  later.  If  boundary  layer  interference  occurs,  the  convective  heat  transfer 
coefficient  between  the  cylinder  and  the  fluid  may  be  adversely  affected.  The  heat 
transfer  from  a  bank  of  cylinders  may  therefore  deviate  from  the  heat  transfer  from  an 
isolated  cylinder  thus  resulting  in  reduced  performance  effectiveness. 
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Figure  4:  Boundary  Layers  on  Cylinders  in  Cross  Flow 


Depending  on  design  constraints,  it  may  become  necessary  to  place  the  tubes  of  a 
heat  exchanger  close  enough  to  each  other  so  that  this  boundary  layer  interference  could 
become  an  issue.  If  this  is  the  case,  the  quantitative  influence  of  boundary  layer 
interaction  on  the  convective  heat  transfer  coefficient  would  need  to  be  known  in  order  to 
more  effectively  design  a  heat  exchanger  for  a  thermoacoustic  refrigerator.  Finding  this 
affect  of  boundary  layer  interference  on  the  convective  heat  transfer  coefficient  is  the 
ultimate  goal  of  this  experimental  study. 


III.  EXPERIMENT 


A.  INTENT 

Since  a  bank  of  cylindrical  tubes  may  have  a  boundary  layer  interference  effect  on 
it’s  convective  heat  transfer  coefficient,  direct  extrapolation  of  experimental  results  for  an 
isolated  cylinder  may  not  be  realistic.  This  experiment  is  designed  to  examine  the  effects 
of  spacing  between  the  cylinders  of  a  tube  bank.  A  long  cylindrical  sound  chamber  was 
chosen  for  this  experiment  because  it  helps  to  restrict  the  sound  waves  to  only  the 
longitudinal  direction  of  motion,  thus  preventing  the  occurrence  and  interference  of 
transverse  waves.  Additionally,  the  sound  chamber  resembles  the  structure  of  a 
thermoacoustic  engine,  which  is  the  practical  application  for  which  this  experiment  is 
intended. 

B.  DIMENSIONLESS  PARAMETERS 

In  order  to  link  one  experimenter’s  results  to  another’s,  and  to  be  able  to  apply 
results  to  a  wide  variety  of  situations  other  than  the  experiment  itself,  dimensionless 
parameters  are  necessary.  Common  convective  heat  transfer  correlations  for  uniform  flow 
over  bluff  bodies  have  long  used  such  parameters  as  the  Reynolds,  Nusselt,  Prandtle,  and 
Grashof  numbers.  Oscillating  flow  is  more  complicated,  however,  and  additional 
dimensionless  parameters  are  necessary  to  fully  describe  the  flow.  Richardson  (1967) 
gives  a  detailed  derivation  of  these  parameters. 
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1.  Cylinder  Length  Scale 

In  order  to  assume  that  the  flow  around  the  tube  bank  is  incompressible,  the 
radian  wavelength  of  the  acoustic  field  must  be  large  compared  to  the  characteristic 
length  of  the  cylinder.  This  characteristic  length  is  chosen  as  the  cylinder’s  diameter,  d, 
and  the  radian  wavelength,  Xr,  is  defined  as 


c 

co 


(1) 


The  requirement  is  that 


(2) 


So  %  is  defined  as  the  ratio  of  the  characteristic  length  of  the  cylinder  to  the  radian 
wavelength 


* 


dco 

= —  « 1 
c 


(3) 


By  ensuring  this  requirement,  Lighthill  (1963)  showed  that  radiation  affects  due  to 
acoustic  streaming  are  very  small  and  can  be  neglected,  making  only  the  acoustic  field 
important. 
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2.  Amplitude  Parameter 


Another  dimensionless  parameter  of  importance  is  s,  defined  as  the  ratio  of  the 
displacement  amplitude  of  particle  oscillation  to  the  characteristic  body  length.  After 
manipulation,  s  can  be  expressed  in  terms  of  the  pressure  ratio  of  oscillation  as  shown  by 
Gopinath  and  Harder  (2000): 


s  = 


c 

doxy 


V  P*  ) 


(4) 


If  s  is  small  compared  to  unity,  the  flow  will  remain  attached  to  the  cylinder  and  will  be 
laminar.  If  s  becomes  greater,  the  flow  will  shed  from  the  cylinder. 

3.  Frequency  Parameter 

A  frequency  parameter,  A  ,  is  defined  as  follows: 


v 


(5) 


Gopinath  and  Harder  (1999)  showed  that  if  A2  is  much  greater  than  unity,  the  Stokes 
shear  layer  is  confined  to  a  narrow  region  and  the  acoustic  streaming  effect  appears  as 
slip  velocity  along  the  cylinder  surface. 
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4. 


Interference  Parameter 


For  cylinders  arranged  with  axes  parallel  and  with  uniform  separation,  let  St  be 
the  center  to  center  distance  between  the  cylinders.  Then  the  distance  between  the  outer 
edges  of  the  cylinders  at  their  closest  point,  is  therefore 

g  =  ST-2r  =  ST-d  (6) 

Riley  (1965)  and  Stuart  (1966)  showed  that  the  boundary  layer  thickness,  8,  of  the  steady 
streaming  flow  generated  on  the  cylinder  at  its  thickest  is  given  by: 


5 


V  co 


(7) 


which  is  used  here  as  the  separation  distance  criterion  for  interference  between  the 
boundary  layers.  For  parallel  cylinders,  boundary  layer  interference  will  occur  if 


(8) 


Substituting  the  approximation  for  8  and  squaring  both  sides,  it  follows  that 


400l> 


<1 


(9) 
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is  the  requirement  for  boundary  layer  interference.  This  is  called  the  interference 
parameter,  cp,  and  can  also  be  expressed  in  terms  of  other  parameters 


400u 


400o 


400 


'N2 


zJL  —  \ 
d 


(10) 


Boundary  layer  interference  will  occur  if  <p  <  1 . 

5.  Streaming  Reynolds  Number 

The  Reynolds  number  is  a  familiar  dimensionless  parameter  in  fluid  dynamics 
and  convective  heat  transfer.  In  oscillating  flow,  however,  the  parameter  is  modified  to 
include  the  frequency  and  pressure  ratio.  Stuart  (1966)  and  Riley  (1966)  found  that  both 
the  frequency  and  amplitude  parameter  are  necessary  to  define  a  suitable  Reynolds 
number.  The  streaming  Reynolds  number,  Rs,  is  defined  as 


Rs=£2  A2=. 


covy 


(11) 


For  buoyancy  effects  to  be  negligible,  Gopinath  and  Harder  (1999)  showed 


(12) 


This  criterion  ensures  that  free  convection  is  small  compared  to  forced  convection. 
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c. 


APPARATUS 


The  equipment  utilized  in  this  experiment  consists  of  three  main  assemblies:  the 
sound  chamber,  the  tube  bank,  and  the  instruments  used  for  data  collection. 

1.  Sound  Chamber 

A  long  cylindrical  tube  is  used  to  obtain  a  standing,  resonant,  isolated  acoustic 
wave.  The  chamber  is  made  of  transparent  plexi-glass.  The  inside  diameter  is  3” 
(~76mm),  and  the  wall  thickness  is  V"  (~6mm).  The  length  of  the  chamber  is  2m.  The 
Chamber  is  supported  by  eight  stantions  along  its  length,  connected  to  a  plexi-glass  base 
plate.  At  one  end  of  the  chamber  is  the  acoustic  driver,  and  at  the  other  end  is  an 
adjustable  endplate  that  allows  the  length  of  the  tube  to  be  varied.  The  end  plate  makes  an 
airtight  seal  to  the  inside  chamber  wall  with  a  rubber  O-ring.  Figure  5  shows  the  sound 
chamber.  The  acoustic  driver  is  a  JBL®  Model  2490H  midrange  compression  driver  with 
8Q  impedance  and  a  frequency  range  of  250Hz  to  2.5kHz. 


Figure  5:  Experimental  Sound  Chamber 
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2. 


Tube  Bank 


The  tube  bank  consists  of  a  central  heated  cylinder  and  two  “dummy”  cylinders 
on  each  side  of  it.  The  heated  cylinder  is  a  smooth  surfaced  Watlow®  cartridge  heater 
with  a  1/8”  diameter.  In  addition  to  the  internal  heating  element,  the  heater  contains  a 
type  J  thermocouple  to  measure  the  surface  temperature  of  the  cylinder.  When  inside  the 
sound  chamber,  the  cylinder  length  to  diameter  ratio  is  ~24,  which  ensures  that  end 
effects  can  be  considered  negligible  for  heat  transfer  calculations.  The  “dummy” 
cylinders  are  smooth  1/8”  diameter  brass  tubes.  The  cylinders  are  put  into  the  sound 
chamber  through  oblong  grooves  cut  through  the  chamber’s  top  and  bottom  surfaces.  The 
grooves  allow  the  spacing  between  the  cylinders,  Sr,  to  be  varied.  The  tubes  are  held  in 
place  by  an  aluminum  plate  fastened  to  the  outside  of  the  chamber.  Figure  6  shows  the 
tube  bank  arrangement.  In  order  to  ensure  that  the  sound  chamber  is  sealed  and  isolated 
from  the  outside  air,  Silly  Putty®  is  used  to  fill  the  grooves  around  the  cylinders. 


Ftrttner 
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3. 


Electronic  Instruments  Package 


A  Hewlett  Packard®  33 120 A  variable  waveform  generator  is  used  to  power  the 
acoustic  driver  with  a  sinusoidal  signal  at  the  desired  frequency.  Prior  to  being  input  to 
the  acoustic  driver,  the  signal  is  amplified  by  a  Techron®  7540  power  amplifier. 

The  adjustable  endplate  of  the  sound  chamber  was  fitted  with  an  Endevco®  model 
8510B-5  pressure  transducer.  The  transducer’s  sensitivity,  S,  is  50.89mV/psi.  The 
transducer  output  goes  to  a  pre-amp  with  a  gain,  G,  of  100,  which  gives  an  amplified  gain 
of  5089m V/psi.  The  pre-amp  then  is  input  into  a  Hewlett  Packard®  3562A  dynamic 
signal  analyzer.  The  signal  analyzer  allows  the  output  from  the  pressure  transducer,  VmjC, 
to  be  observed  in  both  time  domain  and  power  spectrum  outputs.  Using  the  power 
spectrum  output,  the  pressure  ratio  of  the  resonant  wave  inside  the  chamber  can  be 
calculated  as  follows 


(13) 


Where  Pm  is  the  mean  atmospheric  pressure  of  14.7psi,  and  P0  is  the  difference  between 
the  maximum  pressure  in  the  sound  chamber  (which  occurs  at  pressure  antinodes,  or 
velocity  nodes)  and  Pm.  The  power  spectrum  also  allows  disturbances  from  harmonics 
and  other  sources  to  be  observed  so  that  they  can  be  minimized. 

The  power  to  the  heated  cylinder  is  supplied  by  a  Kikusui®  model  PAR  160 A 
regulated  DC  power  supply.  To  measure  the  power  to  the  heater,  the  voltage  drop  across 
it  and  the  current  running  through  it  are  measured.  A  Hewlett  Packard®  34401 A  digital 
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multimeter  is  placed  in  parallel  with  the  heater  to  measure  the  voltage  across  it,  Vh-  The 
current  through  the  heater,  Ih,  is  calculated  by  measuring  the  voltage  drop  across  a 
resistor  of  known  resistance,  Rr,  and  applying  Ohm’s  law 


/  =/ 

H  R  R» 


(14) 


Figure  7  is  the  electrical  circuit  used  to  calculate  the  power  to  the  heated  cylinder. 


Figure  7:  Heater  Circuit 

Finally,  a  Keithley®  740  system-scanning  thermometer  is  used  to  measure  the 
heated  cylinder  surface  temperature  and  the  ambient  fluid  temperature.  As  stated  earlier, 
the  heater  is  supplied  with  an  internal  type  J  thermocouple.  The  ambient  temperature  is 
taken  with  a  type  E  thermocouple  inside  a  small  probe.  The  probe  in  inserted  into  the 
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sound  chamber  via  a  removable  access  port  in  the  side  of  the  chamber.  The  thermometer 
can  be  toggled  between  the  thermocouples  to  record  both  temperatures. 

Figure  8  shows  the  schematic  of  the  instruments  used  in  the  experiment.  Figure  9 
is  a  photograph  of  the  experimental  apparatus. 


Heater 


Figure  8:  Instruments  Schematic 


Figure  9:  Photograph  of  Experimental  Apparatus 


18 


D.  PROCEDURE 


Before  data  “runs”  could  begin,  frequencies  and  corresponding  chamber  lengths 
had  to  be  found  that  would  set  up  a  resonant,  standing  wave  inside  the  chamber  with  a 
velocity  antinode  at  the  tube  bank  location.  From  Figure  10,  it  is  obvious  that  the  distance 


L  must  be  an  odd  multiple  of  A./4.  Put  another  way. 


c 


(15) 


and  n  must  be  an  odd  integer.  To  accomplish  this,  a  pressure  transducer  was  placed  at  the 
chamber  surface  at  the  location  of  the  tube  bank.  The  acoustic  driver  was  then  powered 
with  an  arbitrary  frequency,  and  an  arbitrary  tube  length.  The  frequency  was  adjusted 
until  Vmic  was  at  a  minimum.  Next,  n  was  calculated  to  see  if  it  was  indeed  an  odd 
integer.  If  not,  the  chamber  length  was  adjusted  by  moving  the  endplate  and  the  process 
was  repeated.  In  all,  6  resonant  frequencies  were  found.  These  frequencies  are  listed  in 
Table  1. 


f  (Hz) 

L(m) 

n 

maxPo/Pm  (%) 

383 

0.24 

1.069 

~3.3 

676 

0.65 

5.109 

-3.1 

723 

0.6 

5.044 

-3 

925 

0.65 

6.991 

-2.8 

1202 

0.65 

9.085 

-2 

1461 

0.65 

11.04 

-1.3 

Table  1:  Resonate  Frequencies  Used 
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Figure  10:  Geometry  for  Maximum  Velocity  at  Tube  Bank 

Once  the  frequencies  are  identified,  it  is  necessary  to  determine  the  maximum 
possible  pressure  ratio  for  each  frequency,  without  too  much  interference  from  harmonics 
and  other  sources.  This  is  done  by  observing  the  power  spectrum.  Since  the  first 
harmonic  of  the  frequency,  2f,  is  out  of  phase  with  the  primary  signal,  it  will  not 
interfere.  However,  the  second  harmonic,  3f,  is  in  phase  and  will  therefore  interfere  with 
the  primary  signal.  The  pressure  ratio  was  increased  until  the  transducer  voltage  from  the 
second  harmonic  was  5%  of  the  primary  signal.  This  value  of  5%  interference  was  taken 
as  the  cutoff  for  the  maximum  pressure  ratio.  The  maximum  pressure  ratio  for  each 
frequency  is  listed  is  Table  1. 

After  the  resonant  frequencies  and  the  maximum  possible  pressure  ratios  were 
found,  the  data  collection  could  begin.  The  first  step  was  to  set  the  tube  bank  spacing. 
Four  values  of  S-r/d  were  used  in  this  experiment;  1.25, 1.5, 1.75,  and  2.0.  The  spacing 
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was  measured  to  .001”  using  calipers.  After  the  spacing  was  set,  the  tubes  were  inserted 
into  the  chamber  and  sealed  with  putty. 

Once  the  spacing  was  set,  the  data  collection  could  begin.  First,  the  power  to  the 
heated  cylinder  was  engaged.  The  heater  would  then  begin  to  rise  in  temperature,  which 
was  monitored  by  the  thermometer.  Next,  the  acoustic  driver  was  powered  by  turning  on 
the  power  amplifier.  The  power  to  the  driver  was  increased  until  the  temperature  of  the 
heater  would  level  off  and  attain  a  steady  state  value.  If  the  pressure  ratio  of  the  driver 
was  too  small,  the  heater  would  continue  to  rise  in  temperature.  If  the  pressure  ratio  was 
too  large,  the  heater  temperature  would  begin  to  decrease.  Once  a  steady  state 
temperature  was  achieved,  the  data  was  recorded.  VmiC  was  recorded  first,  then  Vh  and 
Vr.  The  steady  state  temperature  of  the  heater,  Th,  was  recorded,  then  the  power  to  the 
driver  and  the  power  to  the  heater  were  simultaneously  disengaged.  The  ambient 
thermocouple  probe  was  then  inserted  into  the  chamber  and  the  ambient  temperature 
recorded.  This  completed  the  data  collection  for  one  data  point.  For  each  frequency,  at 
each  tube  spacing,  twenty  data  points  were  taken.  The  voltage  across  the  heater  was  then 
raised  slightly  and  the  process  repeated. 

Once  the  data  was  collected,  the  Nusselt  number  and  Streaming  Reynolds 
Number  were  calculated  using  a  computer  spreadsheet.  The  Nusselt  number,  Nua,  is 
defined  as 


Nud 


(16) 


Where  the  convective  heat  transfer  coefficient,  h,  is  defined  as 
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(17) 


h  = 


Qh 

Ah^T) 


The  power  to  the  heater  is 


< 


Qh=IhVh=^~Vh 

-K/? 


(18) 


Additionally, 


Ah  =  7vdLH 

(19) 

(at)=t„-ta 

(20) 

Using  Eqns.  (16)  through  (20),  Nud  was  calculated  by 


Nuj  = 


YZh 

tzRrLh  (Th  —  Ta  )kajr 


(21) 


The  values  of  s,  %,  A2,  (p,  and  Rs  were  then  calculated  as  described  previously.  The  speed 
of  sound,  the  kinematic  viscosity,  and  kajr  were  calculated  at  the  average  fluid 
temperature 


T  +T 

y  _ 


(22) 
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IV.  RESULTS  AND  DISCUSSION 


A  total  of  over  450  individual  data  points  were  calculated  for  the  six  frequencies 
and  four  tube  spacing  values  as  described  in  the  previous  chapter.  Table  2  shows  the 
ranges  of  the  different  dimensionless  parameters  covered  during  the  experiment. 


Parameter 

min 

max 

epsilon 

warn 

i 

X 

0.02 

0.08 

AA2 

1500 

5900 

phi 

0.2 

15 

Po/Pm 

0.50% 

3.30% 

Rs 

<10 

1500 

Table  2:  Parameter  Ranges  Used 


For  values  of  (p>l,  the  results  replicate  those  of  an  isolated  cylinder  (to  within 
experimental  error).  Although  this  conclusion  might  have  been  expected  since  there  is  no 
boundary  layer  interference  anticipated  in  this  regime,  it  nonetheless  provides  a  firm 
corroboration  with  the  earlier  single  cylinder  results  of  Gopinath  &  Harder  (2000).  For 
“low”  values  of  Rs  the  relationship  between  Nua  and  Rs  is  expected  to  be  of  the  form 


Nud  =  FPr*  Rsos 


(23) 


which  exemplifies  the  laminar,  attached  flow  regime.  This  low-amplitude  attached  and 
laminar  behavior  of  the  flow  has  also  been  demonstrated  by  Sarpkaya  (1986)  in  an 
independent  context,  i.e.  on  the  basis  of  force  measurements  on  a  cylinder  in  an 
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oscillatory  flow.  Since  Pr  is  constant  throughout  this  heat  transfer  experiment  (Prair=0.7), 
it  can  be  included  in  the  constant  Y  and  the  relationship  becomes 

Nud=CR° 5  (24-a) 

Gopinath  and  Harder  (2000)  found  that  at  larger  values  of  Rs,  the  data  has  a 
steeper  trend,  roughly  proportional  to  Rs° 75  which  may  be  attributable  to  vortex  shedding 
resulting  from  flow  instabilities.  For  Rs  values  >500,  the  data  was  correlated  using  the 
form 


Nud  =  CRS015  (24-b) 

Table  3  shows  the  values  of  the  constant  C  in  Eqns.  (24-a)  and  (24-b),  correlated  using  a 
transition  of  Rs=500.  Figures  1 1-13  are  logarithmic  plots  of  Nuj  vs.  Rs  for  Sx/d  values  of 
1.5, 1.75,  and  2.0,  respectively.  Best-fit  lines  are  constructed  from  Eqns.  (24-a)  and  (24- 
b).  Error  bars  from  an  uncertainty  analysis  (Appendix  C)  are  included  for  arbitrary  data 
points.  Figure  14  is  a  plot  of  Nuj  vs.  Rs  for  all  data  corresponding  to  cp>l. 

Davidson  (1973)  predicted  a  correlation  of  Nud=1.05Rs° 5  for  the  attached  flow 
regime,  while  Gopinath  and  Harder  (2000)  found  the  leading  coefficient  to  be  0.90.  For 
“large”  values  of  Rs,  Gopinath  and  Harder  (2000)  have  suggested  that  the  correlation 
should  be  Nud=0.20Rs°'75.  Table  3  also  includes  the  percentage  deviation  of  the  current 
experiment  from  these  previous  correlations. 
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Rs<500 


Rs>500 


Figure  11:  Nusseit  Number  vs.  Streaming  Reynolds  Number  for  ST/d=1.5. 

The  solid  line  is  a  Rs0'5  fit  through  the  data.  The  dashed  line  is  a  Rs  0/75  fit  in  the  vortex  shedding 
regime.  Error  bars  are  included  for  arbitrary  data  points. 


Figure  13:  Nusselt  Number  vs.  Streaming  Reynolds  Number  for  ST/d=2.0. 

The  solid  line  is  a  Rs°‘5  fit  through  the  data.  The  dashed  line  is  a  Rs  075  fit  in  the  vortex  shedding 
regime.  Error  bars  are  included  for  arbitrary  data  points. 


Figures  11-14  indicate  that  the  data  from  individual  frequencies  cannot  be 
distinguished  from  one  other.  With  cp>l,  the  data  falls  along  roughly  the  same  curve.  Nud 
is  independent  of  the  tube  spacing,  and  <p  is  not  necessary  to  correlate  the  data. 

Figure  15  is  a  logarithmic  plot  of  Nud  vs.  Rs  for  Sx/d=1.25.  For  this  close  spacing, 
(p<l  for  all  frequencies  used,  and  boundary  layer  interference  becomes  an  issue.  From 
Figure  15  it  is  obvious  that  the  data  for  the  different  frequencies  do  not  fall  on  the  same 
curve.  Nud  can  not  be  correlated  with  Rs  alone,  and  the  influence  of  cp  must  also  be 
included.  This  data  was  correlated  using  the  form 

Nud  =  C<p°Rsh  (25) 

Again,  there  was  a  transition  at  Rs=500,  so  the  correlations  used  b=0.5  for  Rs<500,  and 
b=0.75  for  Rs>500.  The  following  results  were  found 


Nud  =  l.07(p019  Rs0'5 

Rs<500,  cp<l 

(26-a) 

Nud  =  0.21<pou  Rs°'15 

Rs>500,  cp<l 

(26-b) 

Eqns.  (26-a)  and  (26-b)  show  that  as  cp-»l,  the  correlations  converge  to  those  of  an 
isolated  cylinder. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Experiments  were  performed  to  determine  the  effect  of  tube  spacing  on 
convective  heat  transfer  from  a  tube  bank  in  a  zero  mean  oscillating  flow.  A  tube  bank  is 
used  to  model  a  heat  exchanger.  The  role  of  the  interference  of  boundary  layers  on 
adjacent  tubes  in  the  tube  bank  was  discussed.  The  effect  of  boundary  layer  interference 
is  needed  to  design  more  efficient  heat  exchangers  for  thermoacoustic  engines. 

Without  boundary  layer  interference,  the  heat  transfer  from  a  tube  bank  closely 
follows  that  of  an  isolated  cylinder.  It  was  found  that,  as  expected,  the  Nusselt  number 
does  not  depend  on  tube  pitch  for  sufficiently  large  spacings  corresponding  to  values  of  9 
greater  than  1.  However  for  values  of  9  less  than  1,  the  Nusselt  number  must  be 
correlated  to  9  as  well  as  Rs.  As  the  tube  spacing  and  hence  9  decreases,  the  heat  transfer 
rate  characterized  by  Nud  is  degraded.  However,  the  decrease  in  Nua  is  not  great, 
especially  at  high  values  of  Rs.  At  high  streaming  Reynolds  numbers,  9  can  be  as  low  as 
0.12  although  Nud  is  still  at  80%  of  its  value  for  9=1 .  At  low  streaming  Reynolds 
numbers,  9  can  be  as  low  as  0.32  and  still  have  less  than  20%  degradation  in  Nud.  It  can 
be  concluded  that  only  a  small  sacrifice  in  heat  transfer  will  result  from  having  closely 
spaced  tubes.  Provided  that  9  is  not  too  small,  heat  exchanger  designs  can  use  a  cluster  of 
many  closely  spaced  tubes,  allowing  more  cooling  fluid  to  be  transported,  without 
significantly  degrading  the  heat  transfer  characteristics  of  each  tube. 

Based  on  a  random  error  propagation  analysis,  it  was  determined  that  the  greatest 
source  of  uncertainty  came  from  temperature  measurements.  Future  experiments  would 
benefit  from  more  precise  measurements  of  temperature.  Additionally,  the  way  in  which 
the  ambient  temperature  is  measured  could  be  improved.  Instead  of  inserting  a 
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thermocouple  probe  into  the  sound  chamber  after  the  data  run,  a  thermocouple 
permanently  fixed  to  the  inside  surface  of  the  sound  chamber,  in  such  a  way  as  to  not 
interfere  with  the  flow  characteristics,  would  allow  the  ambient  temperature  and  cylinder 
temperature  to  be  recorded  simultaneously. 

To  carry  this  study  further,  it  is  recommended  that  a  similar  experiment  be 
conducted  to  determine  the  affect  of  an  in-line  arrangement  of  cylinders,  that  is,  a  bank  of 
cylinders  arranged  with  their  plane  parallel  to  the  direction  of  fluid  oscillation.  Such  an 
experiment  could  study  both  aligned  and  staggered  cylinders  as  in  conventional  heat 
exchanger  tube  bank  arrangements.  Additionally,  determining  the  heat  transfer 
characteristics  of  finned  cylinders  would  greatly  benefit  heat  exchanger  design.  Finally, 
experiments  with  fluids  other  than  air  would  allow  a  Prandtl  number  dependence  to  be 
studied. 
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APPENDIX 


A.  SAMPLE  CALCULATIONS 

The  following  sample  calculations  use  the  following  data: 
Measured  Values: 


ST/d=1.25 
L=65  m 
f=676  Hz 
Vmic=1878  mV 
VR=0. 17225  V 
VH=1 0.2805  V 
Th=28.6  °C 
Ta=22.9  °C 


Constants: 


Rair=287  m2/s2K 
7=1.4 

Pm=14.7psi 
Ln=.076m 
S=50.89  mV/psi 
G=100 
d=.003175m 


1.  Nusselt  Number 

The  average  fluid  temperature  is: 


T  = 

avg 


Tu+Ta  _28.6°C  +  22.9°C 


25.75°C  =  298.9K 


From  this,  the  following  values  are  taken  from  Kays  (XXX): 
kair=.0260646  W/mK 
v=.00001563  m2/s 
AT  is  calculated  as  follows: 


(A- 
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The  Nusselt  number  is  then  calculated: 


Nu,= 


vrvh 


(0.17225FXiO-2805F) 


7rRRLH{hJ)kair  ^(2QX-076wX5.7^)(.0260646^/  J 


=  24.96  (A-3) 


2.  Streaming  Reynolds  Number 

The  speed  of  soimd  is  given  by: 

C  =  =  346.55%  (A-4) 


The  pressure  ratio  is  given  by: 


1878 mV, 


K.„ 


4SG 


50.89wF/ 


>)10°) 


14.7/w/ 


=  .0251 


(A-5) 


The  angular  frequency  is: 

a  =  2nf  =  [in  ra%,c/e)(676flz) = 4247.433 rad/s  (A-6) 


The  streaming  Reynolds  number  is  then: 


RS  = 


c2 

f  P.) 

2 

=  l 

1 

(346-55%') 

f 

_ i _ 

cooy1 

(4247.433™%’ 

(.00001 563 mYs) 

ti-4 y 

(.0251)2  =581.5  (A-7) 


B.  UNCERTAINTY  ANALYSIS 


This  uncertainty  analysis  is  done  using  uncertainty  propagation  found  in  Beckwith 
etal.  (1993). 

1.  Nusselt  Number 

The  Nusselt  number  equation  is: 


Ar  hd  VHVR 
Nu,  =  ■  -  H  R 


Kir  * LHRRATkair 


(A-8) 


The  fractional  imcertainty  of  the  Nusselt  number  is: 


u 


Nu, 


Nud  \ 


r 

AT. 


uv 


V 


yVRj 


uv 


V  ru.  V 


\VH  J 


\rr  j 


Ur 

\Lh  j 


V  (u„  \ 


_Kir 

V  Kir  J 


(A-9) 


The  uncertainties  of  the  temperature  measurements  are  0.5°C.  With  AT=Th-Ta, 
the  imcertainty  of  the  temperature  difference  is: 


UA t  ~  ■ 


SAT 

KdTH 


^2  ( 


■Ur 


dAT 

- UT 

dTA  T*j 


=  V[(lX0-5°C)f +  [(lX0.5°C)]2  =  .7071  °C  =  .7071^ 


(A- 10) 


The  uncertainties  in  the  voltage  measurements  are  taken  from  the  HP 
multimeter’s  users  manual: 

uvR=uvH=0.05%reading+0.02%full  scale  output 
The  resistance  uncertainty  is  1%  of  the  resistor  value.  The  uncertainty  of  the  heater  length 
is  0.001”.  The  uncertainty  of  kairis  0.00005  W/mK. 

2.  Streaming  Reynolds  Number 

The  streaming  Reynolds  number  is: 
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Rs  = 


R 


T  V  . 

avg  nuc 


K2nypm  j 

The  fractional  uncertainty  is  then: 


2  o2 


fuGzS 


(A-ll) 


Ur 


—  = 


l'UT..)2  i'"/  Y  fu..Y  1 


T 

V  lovz  J 


f 


\0  J 


Uy 

2_jw_ 

v  K,ic  j 


+[  2— 
l  G 


v  r  uy 

2— 

)  V  S  J 


(A- 12) 


With  Tavg=(TH+TA)/2,  the  uncertainty  in  Tavg  is: 


ur  = . 

lavg 


dT 


avg 


Ur 


V  dTH  J 


dT 


V 


avg 


ydTA 


f 0.5°C^2  f 


J 


0.5  °C 


=  0.3536°C  =  0.3536AT 


(A- 13) 


The  uncertainty  in  the  frequency  is  taken  from  the  HP  waveform  generator  user’s 
manual: 

uj=20x10'6x  reading 

The  uncertainty  in  u  is  5x1  O'10  m2/s.  The  uncertainty  in  VmjC  is  5mV.  The  uncertainty  in  G 
is  0.4.  The  uncertainty  in  S  due  to  nonlinearity  is  given  by  the  pressure  transducer 
calibration  data  as  0.381mV/psi. 
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c. 


EXPERIMENTAL  DATA 


The  following  pages  contain  the  data  gathered  through  experimentation  in 
computer  spreadsheet  form.  All  important  parameters  are  listed,  although  some  constants 
have  been  omitted  to  reduce  size. 
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UdT/T  lUVR/VR  lUVHA/H  |Uk/k  |UNu/Nu  |UNu  I  lUTavg/Tavg  |Uf/f  lUv/v  |Uvmic/Vmic  IUG/G  I  US/S  |URs/Rs  [URs 


39 


f  (Hz)  IVmic  (mV)  (VR  (V)  |VH  (V)  |TH  (C)  |TA(C)  IdeltaT(C)  Itavg  (C)  Jk(W/mK)  INu  |c(m/s)  |v  (m*2/s)  I  PR  (epsilon  |X  [7^2 


If  (Hz)  |Vmic(mV) |VR(V)  |VH  (V)  |TH  (C)  |TA(C)  jdeltaT(C)  |Tavg(C)  Jk(W/mK)  |Nu  |c(m/s)  lv(mA2/s)  (PR  (epsilon 


UdT/T  l  UVR/VR  lUVH/VH  |Uk/k  |UNu/Nu  |UNu  1  lUTavgfTavg  |Uf/f  |Uv/v  |Uvmic/Vmic  lUG/G  I  US/S  |URs/Rs  |URs 


41 


42 


if  (Hz)  (Vmic(mV)  |VR(Vj  |VH(Vj  |TH(C)  |TA(C)  IdeltaT(C)  ITavg(C)  |k(W/mK)  |Nu  |c(m/s)  |v(m*2/s)  |PR  lepsiion  IX  p2 


UdT/T  lUVR/VR  lUVHA/H  |Uk/k  |UNu/Nu  |UNu  I  lUTavgrfavg  |Uf/f  |Uv/v  |UvmicA/mic  lUG/G  I  US/S  lURs/Rs  lURs 


43 


44 


Ivmic(mV)  VR  (V)  |VH  (V)  |TH(C)  |TA(C)  IdeltaT(C)  iTavg  (C)  |k(W/mK)  |Nu  |c(m/s)  |v{mA2/s)  |PR  I  epsilon 


UdT/T  IUVRA/R  lUVHA/H  |Uk/k  |UNu/Nu  |UNu  I  iUTavgnavg  |Uf/f  |Uv/v  |Uvmic/Vmic  lUG/G  IUS/S  |URs/Rs  |URs 
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46 


UdT/T  lUVR/VR  IUVH/VH  |Uk/k  |UNu/Nu  |UNu  |  |UTavg/Tavg  |Uf/f  |Uv/v  |Uvmic/Vmic  IUG/G  IUS/S  lURs/Rs  |URs 


47 


48 


If  (Hz)  IVmic(mV)  |VR  (V)  |VH  (V)  |TH  (C)  |TA(C)  IdeltaT(C)  |Tavg(C)  |k(W/mK)  |Nu  |c(m/s)  |v{mA2/s)  |PR  lepsilon 


UdT/T  IUVRA/R  lUVH/VH  IUkA<  |UNu/Nu  |UNu  |UTavg/Tavg  |Uf/f  [Uv/v |UvmicA/mic  IUG/G  lUSTs  |URs/Rs  lURs 


49 


UdTYT  IUVRA/R  lUVH/VH |Uk/k |UNu/Nu  |UNu  I  |UTavg/Tavg  |Uf/f  |Uv/v  |UvmicA/mic  |UG/G  | US/S  |URs/Rs  |URs 


51 


52 


UdT/T  1UVR/VR  lUVHA/H  |Uk/k  lUNu/Nu  |UNu  I  |UTavg/Tavg  |Uf/f  |Uv/v  |UvmicA/mic  lUG/G  lUS/S  |URs/Rs  |URs 
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54 


If  (Hz)  IVmic(mV) [VR(Vj |VH (V)  |TH  (C)  |TA(C)  IdetlaT(C)  |Tavg  (C)  |k(W/mK)  |Nu  |c(m/s)  lv(mA2/s)  1PR  [epsilon 


0.129  0.00027  0.000121  4E-04 1  0.129  1.26  I  0.02476711  0.000021  3.1E-05I  0.00582751  0.00341  0.007491  0.03193261  3.59951 
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If  (Hz)  |Vmlc(mV)  |VR  (V)  |VH  (V)  |TH  (C)  |TA(C)  IdeUaT(C)  iTavg  (C)  |k(W/mK)  INu |c  (m/s)  |v  (mA2/s)  |PR  lepsilon 


UdT/T  IUVR/VR  lUVH/VH lUk/k  |UNu/Nu  |UNu  I  |UTavg/Tavg  |Uf/f  |Uv/v  |UvmicA/mic  lUG/G  lUS/S  lURs/Rs  |URs 
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58 


UdTTT  IUVR/VR  IUVHA/H  |Uk/k  |UNu/Nu  |UNu  I  lUTavg/Tavg  |Uf/f  |UWv  |UvmicA/mic  IUG/G  I  US/S  |URs/Rs  |URs 


59 


UdT/T  lUVR/VR  lUVHA/H  |Uk /k  |UNu/Nu  |UNu  I  |UTavg/Tavg  |Uf/f  |U vN  lUvmicA/mic  lUG/G  lUsTs  |URs/Rs  lURs 


61 


UdT/T  IUVR/VR  IUVHA/H  |Uk/k  lUNu/Nu  |UNu  |  lUTavg/Tavg  |Uf/f  [Uv/v  |UvmicA/mic  IUG/G  I  US/S  |URs/Rs  |URs 


63 


0.1607  0.00027  0.0001 1  0.00038 1  0.161021  1.90091  I  0.024810531  0.00002  3.  IE-06  0.0049116  0.0034  0.0074867  0.03134 
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68 


If  (Hz)  Ivmic(mV)  IvR(V)  IVH  (V)  |TH(C)  JTA(C)  IdeltaT(C)  ITavg  (C)  |k(W/mK)  INu  jc(m/s)  |v  (mA2/s)  IPR  lepsiton  IX  p2 


UdT/T  lUVR/VR  IUVHA/H  |Uk/k  lUNu/Nu  |UNu  I  |UTavg/Tavg  |Uf/f  lUv/v  |UvmicA/mic  lUG/G  I  US/S  |URs/Rs  |URs 


69 


70 


12021  420  0.069881  4.17140  29.8  26  3.8  27.9  0.026222  6.126]  347.7961  1.58E-05I  0.00561 1  0.058171  0.068951  4810.051  16.27  ]  6.764141 
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